Non-Fermi liquid behavior due to U(l) gauge field in two dimensions 
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We study the damping rate of massless Dirac fermions due to the U(l) gauge field in (2+1)- 
dimensional quantum electrodynamics. In the absence of a Maxwell term for the gauge field, the 
fermion damping rate ImS(w, T) is found to diverge in both perturbative and self-consistent results. 
In the presence of a Maxwell term, there is still divergence in the perturbative results for Im£(a>, T). 
Once the Maxwell term is included into the self-consistent equations for fermion self-energy and 
vacuum polarization functions, the fermion damping rate is free of divergence and exhibits non- 
Fermi liquid behavior: ImE(w, T) oc max( v / aJ, y/T). 

PACS numbers: 



I. INTRODUCTION 

In a normal Fermi liquid, the Landau quasiparticles are well-defined in the low-energy regime since their damping 
rate vanishes rapidly as ImE(w,T) oc max(oj 2 ,T 2 ) upon approaching the Fermi surface. Such rapidly diminishing 
damping rate is guaranteed by the Pauli exclusion principle and can be naturally understood by the fact that most of 
the states into which the fermions on Fermi surface would be scattered are already occupied by other fermions. Indeed, 
the states on the Fermi surface have an infinite lifetime and the low-lying quasiparticles can be safely considered as 
being nearly independent even when the interaction is not weak. Generally, one condition for the stability of Fermi 
liquid is the absence of singular or long-range interaction [l[. In normal metals, although the bare Coulomb interaction 
is long-ranged, it becomes short-ranged after including the dynamical screening effect. Therefore, the Fermi liquid 
theory can provide an excellent description for the electron liquid in metals. 

Unlike Coulomb interaction, the long-range gauge interaction usually can not be fully screened in the absence of 
gauge symmetry breaking (it becomes short-ranged in a superconductor via Anderson-Higgs mechanism). It is thus 
possible that the unscreened gauge interaction would generate behaviors those are beyond the scope of Fermi liquid 
theory. In the context of ordinary metals, Holstein et al. first pointed out that non-Fermi liquid behavior arises from 
the coupling of electrons with the unscreened electromagnetic field ■ This prominent result began the adventure for 
seeking various non-Fermi liquid behaviors in several different gauge theories (3rJl6j in the subsequent three decades. 

The (2+l)-dimensional quantum electrodynamics (QELV) of massless Dirac fermions is an interesting model that 
has been widely studied in both high energy physics [T3, El] and condensed matter physics • The gauge field 

is strongly interacting with Dirac fermions, giving rise to rather unusual behaviors. As shown by Appclquist et al., 
the massless Dirac fermions can acquire a finite mass via the mechanism of dynamical chiral symmetry breaking [l7j . 
In the context of condensed matter physics, this phenomenon is usually identified as the formation of long-range 
antiferromagnetism in two-dimensional quantum Hcisenberg antiferromagnet Q. In the phase with chiral symmetry 
unbroken, the Dirac fermions are massless and the gauge field is strongly fluctuating but stays in the deconfmcd 
phase . Now this field theory can describe the physics of a U(l) spin liquid, which is a novel state of matter 
without any spontaneous symmetry breaking. It is important to emphasize that the absence of symmetry breaking 
does not mean the absence of interesting physics. In fact, the U(l) spin liquid manifests a great deal of unusual 
physical properties and has been used to understand several correlated electron systems, including high temperature 
copper-oxide superconductors and some spin- 1/2 Kagome systems [llj]. One of the most remarkable features of 
the gauge interaction is its ability to produce non- Fermi liquid behaviors. Owing to the chirality and linear spectrum 
of Dirac fermions, this interacting field theory is expected to display distinct behaviors compared with the much 
studied non-relativistic gauge systems [l3l - [l6j . 

In this paper, wc study the possible non-Fermi liquid behavior by computing the damping rate ImS(oj, T) of massless 
Dirac fermion within the QED3 theory. Generically, there are two kinds of QED3 theories, depending on whether 
the theory has an explicit Maxwell term or not. First of all, as a well-defined quantum gauge field theory, it contains 
the Maxwell term oc F^F^ explicitly in the action. If the theory is constructed by considering the quantum phase 
fluctuations in underdoped high temperature superconductors, then there is also an explicit Maxwell term in the 
Lagrangianfiol [Tl| . On the contrary, when the effective QED3 theory is obtained by the slave-particle treatment of 
t-J model, there is no Maxwell term in the Lagrangian and the gauge field has its own dynamics only after integrating 
out the matter fields 0-0|- Here, we consider both of these two kinds QED3 and show that they behave differently. 

We first consider the QED3 theory without Maxwell term. In the Coulomb gauge, the temporal and spatial 
components of U(l) gauge field are decoupled. We calculate the longitudinal and transverse fermion damping rates 
at both zero temperature, T = 0, and finite temperature. By straightforward perturbation computation, we find 



2 



that the damping rate is always divergent, cither in the longitudinal contribution or in the transverse contribution. 
The appearance of divergences indicates the insufficiency of ordinary perturbative expansion in treating systems with 
singular interaction. Moreover, divergence still exists in the self-consistcntly coupled equations for fermion damping 
rate and gauge boson propagator. 

We then consider the case with explicit Maxwell term for the gauge field. At the perturbative level, the fermion 
self-energy function diverges even in the presence of such term. After studying the self-consistent equations for fermion 
damping rate and gauge boson propagator, we found that the fermion damping rate is free of divergence and given 
by ImE(o;,r) cx maxfyti, vT), which is a non- Fermi liquid behavior. 

We also discuss the damping rate of Dirac fermions due to the long-range Coulomb interaction. Although Coulomb 
interaction may be identified as the non-relativistic counterpart of the U(l) gauge interaction, they lead to differ- 
ent properties of Dirac fermions. Marginal Fermi liquid behavior is found by both perturbative and self-consistent 
approaches. 

The damping rate of massless Dirac fermion is studied in QED3 without Maxwell term in Sec. 2 and with Maxwell 
term in Sec. 3. In Sec. 4, we consider the Coulomb interaction and discuss why divergence appears in gauge theory, 
but not in such system. We end the paper with a summary and a brief discussion. 



II. FERMION DAMPING RATE IN QED 3 WITHOUT MAXWELL TERM 



The Lagrangian of QED3 has the form 

JV 



C = *I ( d r - iao - iff ■ (8 - *a)) - -F^F^ . (1) 
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In principle, the fermion field can be expressed in four-component or two-component representation. Since we 
only consider the chiral symmetric phase, we adopt the two-component representation of spinor field. The Dirac 
fermion flavor is taken to be a general ./V in order to perform 1/N expansion. In this paper, we discuss only non- 
compact QED3 and hence there are no instantons. Note that the coupling between massless Dirac fermion and 
gauge field respects the Lorentz invariance, so there is no singular velocity renormalization [||. When applied to 
high temperature superconductor, there is indeed an velocity anisotropy. However, this anisotropy turns out to be 
irrelevant, thus restoring the Lorentz invariance [lSj. In this paper, the fermion velocities are simply taken to be 
unity. 

As mentioned in the Introduction, we will first consider QED3 theory without Maxwell term for gauge field. This 
is the effective low-energy theory of t-J model, obtained by using the slave-particle treatment. Now the oc F^ V F^ V 
term is simply dropped from the Lagrangian. The gauge field appearing in such model has its own dynamics only 
after integrating out the fermion fields, as well as other possible matter fields. 

The Matsubara propagator of massless Dirac fermion is 

G (w„,k) = , -, (2) 

tu! n — a ■ k 

where uj n = (2n + 1)ttT with n being integers. After analytic continuation, the retarded propagator reads 

G (w,k) = \— — . (3) 

ui — a ■ k + 10 

For simplicity, the fermion energy is approximated by cr-k ~ |k|. To decouple the temporal and spatial components of 
gauge field, it is convenient to work in the Coulomb gauge /c;<2; = 0. In the imaginary time formalism, the propagator 
for the gauge field can now be written as 

D 00 (fi ro ,q) = 7T77T V ( 4 ) 

^H"m, qj 

where f2 m = IrrmT for bosonic modes with m being integers. The vacuum polarization functions Di(f2 m ,q) and 
Dj_(fi m ,q) come from the one- loop bubble diagram of Dirac fermions to the leading order of 1/N expansion. In 
particular, the polarization function appearing in the spatial component is given by 

D±{Q m ,q) = £> 2 (fi m ,q) - % £>i(fi m ,q). (6) 

q 



3 



The functions _Di(fi m ,q) and £>2(^m,q) arc defined as 

£»i(O ro ,q) = -NTY] [ — 4Tr[GoK,k)Go(O m + W „,q + k)], (7) 

£> 2 (fi mi q) = NTY.J ^^Tr[ ( riGo(a;„,k)cr i Go(O m +w„,q+k)]. (8) 
At zero temperature, T = 0, it is straightforward to show that the retarded polarization functions have the forms 

Di ^ q) = + ^(|n|-|q|) 

16v/|q| 2 - 16V« 2 - |q| 2 

D ± (n,q) = ^0(|q| - |fi|) V|q| 2 - & - i S gnO^0(|f2| - |q|)^ 2 - |q| 2 - (10) 

The fermion damping rate can be calculated by either the Fermi golden rule or the, basically equivalent but more 
formal, diagrammatic many-body technique [2(j. We will utilize the latter one since it is easier to write down the 
self-consistent equations using diagrammatic technique. 

A. Perturbative computation of fermion damping rate 

We now calculate the fermion damping rate using conventional perturbative method. To the lowest order of 1/JV 
expansion, the one- loop self-energy of Dirac fermion is given by Fig. [T] which can be written as 



where 



E(w„,k) = E L (w„,k) + E T (w n ,k), (11) 
S L (w„,k) = -TV / 4-%Go(wn + nm,k + q)D o(nm,q) (12) 

/cPq 
——2<TiG (uj n + O m ,k + q)crj Dij(n m , q) (13) 

which represent the contribution from the longitudinal and transverse gauge field, respectively. The damping rate of 
masslcss Dirac fermion can be obtained by making analytic continuation, iu n — > uj + iS, as 

E(w,k) = E L (w,k) + E T ( tl ;,k), (14) 

and then taking the imaginary part, ImE(o;,k). 



FIG. 1: Fermion self-energy correction to the leading order. The solid line represents the massless Dirac fermion, and the 
wiggly line represents the gauge field. 

We first consider the transverse contribution to the damping rate at T — 0. Using the spectral representations for 
Dirac fermion and gauge boson propagators 

n< , o i , i f + °° duJl ImG (wi,k + q) 

G (uj n + n mi k + q) = - / : — , (15) 



i r +o ° du 2 Im Di(ki) 



D±(tt m7 q) J_ 00 7T iQ m -u] 2 



(16) 
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the imaginary part of retarded self-energy function can be cast in the form 



d 2 q 

(27T): 



+00 



ImExO^k) = 

x [n B (ui — CJ) + n F (bJi)] . 
In the T — limit, the occupation numbers simplify to 



Im[G (wi,k + q)]CTj((Ji. 



9iQj NT 



D±(uii -w,q) 



tib^i — w) + Ui?(wi) — > — — wi). 



(17) 



(18) 



In the above expression, the imaginary retarded fermion propagator is Im[Go(w, k)] = —ttS(uj — |k|), while the photon 
propagator has the form 



Im- 



1 



16 



mi 



|q|). 



D ± (n,q) Nsgnny/Q 2 - |q| : 
After straightforward computation, the transverse fermion damping rate is found to be 



ImEx (uj, k) 



Nir 2 



d 2 q 



|q|)- 



(19) 



(20) 



'( W -|k + q|r-|q|' 

To get an analytic expression, we first choose to use the frequently used on-shell approximation tu = = |k| and get 



ImEx (wk) 



N77 2 



d 2 q 



(|k|-|k + q|) 2 -|q| 2 



(21) 



Since the step function always satisfies 0(|k| — |k + q| — |q|) = 0, we know that ImEx (wk) = 0. More generally, 
6 (u - |k+ q| - |q|) = for |w| < |k|, so that 



ImEx (w< |k|,k) =0. 

In order to get a finite analytic expression, we use the zero-momentum limit |k| = and finally have 

|q| 8^ 



ImEx 0*0 



Nn 



d\q\ 



2w q| 



3iV7T 



(22) 



(23) 



This linear- in-energy expression clearly signals a marginal Fermi liquid behavior 21 1. 

This result is obtained for d > 0. When ui < 0, the expression should be ImEx(w) = jj^u- In general, the 
fermion damping rate has the form ImEx 0*0 = ~~ 3lvxM- Without loss of generality, we consider only positive u in 
the following. The real part of retarded fermion self-energy can be directly obtained by the Kramers-Kronig relation, 
as RcEx(w) awlnw. 

This contribution has its own physical application. In the effective gauge theory derived by the slave-boson approach, 
the gauge field also couples to non-relativistic scalar bosons which describe the motion of charged holons [3, [H-0] . 
The scalar bosons are incompressible and thus effectively screen the temporal component of the gauge field. So the 
temporal component ao can be omitted, but the transverse components a remain unscreened and should be carefully 
treated. Within this effective field theory, several interesting results have been obtained, including the singular 
corrections to specific heat and susceptibility Q, and the algebraic correlation Q. The above results show that the 
transverse gauge field leads to marginal Fermi liquid behavior at T = 0. 

It seems impossible to get an analytical expression for the general damping rate ImEx (w, k), so we compute it by 
numerical skills with the results being shown in Fig. [5J 

The longitudinal contribution to fermion self-energy function will be calculated analogously. The longitudinal 
damping rate can be written as 



ImE L ( W ,k)= / -^SL r^Im[G K,k + q)]Im 
J (2tt) Jo * 



1 



D 1 (wi - w,q) 



Using the expression for temporal gauge propagator 



Im 



1 



Di(fi,q) 



-i6yr> 2 -|q| : 

sgnOTV|q| 2 



■«-|q|), 



(24) 



(25) 





(a) (b) 
FIG. 2: Transverse contribution to damping rate of Dirac fermion without Maxwell term of gauge field at T = 0. 



we obtain 



IiiiEl (w, k) 



Nir 2 



d 2 q 



( W -|k + q|) 2 -|q| 2 



0( w -|k + q|-|q|) 



Analogous to the transverse contribution, we can easily get 

ImS L (lu < |k|,k) = 0. 
In the limit |k| =0, the longitudinal damping rate is 



ImS L (u) 



Ntt 



(26) 



(27) 



(28) 



Obviously there appears a serious infrared divergence. In fact, ImSL (w, k) is always infrared divergent for to > |k| 
The general, ImEL(o;,k) can be written as 



ImEL (lu, k) 



+oo 



rf|q]*i(|q|) 



with 



FiM. 



^ , V^-l k + ql) 2 -lql 2 



dtp 



io |q| 

Where <p is the angle between q and k. When u > |k|, 



lim F x (|q| 

|q|->0 



(w-|k + q|-|q|; 



2tt (w- |k|) 



(29) 



(30) 



(31) 



which is divergent in the infrared region. The general damping rate ImSL(w,k) are shown in Fig. |3l 

Here, an infrared divergence appears in the expression of fermion damping rate due to the singular gauge interaction. 

Such divergent damping rate is surely not well-defined at T — 0. This does not imply that the damping rate itself 

diverges, but reflects the inefficiency of naive perturbation computation. 

We now extend the above consideration to finite temperature, T 3> lu. We first calculate the longitudinal damping 

rate of Dirac fermions. After a series of manipulations, in the limit |k| = ,the longitudinal component of imaginary 

sclf-cncrgy function is written in the form 



ImE L (w,T) = ±- f d|q||q|Imf 1 

SttJo LL>i(|q| -w,q, 



[n B (\q\ -oj)+n F {\q\ 



(32) 



FIG. 3: In the yellow region where u) > |k|, ImTt^iuj, k) is infrared divergent. In the green region where uj < |k|, ImEL(w, k) = 0. 



Notice that the occupation numbers (|q| — ui) and np (|q|) damp exponentially with -^r, so the dominant contri- 



bution of the integral comes from the domain |q| < T. Hence the ultraviolet cutoff of the integral can be set to be T. 
In the high temperature limit, it is convenient to make the simplifications 



n B (|q| - ui) 



T 



|q| -w 



(33) 



The polarization function Z?i(|q| — uj, |q|,T) in the limit T S> u has a very complicated expression, which makes 
analytical computation difficult. To keep the analytic tractability, it is necessary to make some approximations. 
To simplify the analysis, we divide the whole domain of |q| into three sections and then perform the momentum 
integration respectively. Specifically, we decompose the longitudinal damping rate into the following three parts: 



ImS L (cj,T) = — 
= h 




(f|q||q|Im 



1 



Di(|q|-w,q,T) 



[n B {\q\ -u)+n F {\q\ j 



(34) 



Here, the variable q* is a particularly chosen quantity between ui and T which we employ to simply the calculation. 
The key motivation to introduce this quantity is to specify the most important contribution of the integration. In 
principle, we can employ any q* between uj and T . In this paper, we assume that q* satisfy 



(35) 




The polarization functions are given by the following expressions. When < |q| < 

T(w-|q|) 



R&Di(|q|-w,|q|,T)«CiT-aj 



y/u 2 -2u\q\ 



ImDidql-w.lql.T) w-C a 



|q| 2 ( W - 



Ty/U 



2uj q 



(36) 
(37) 



When % < |q| < q* 



Re£>i(|q|-w,|q|,r) 
LnDidql-w.lql.r) 



CiT + C 4 
C 2 



T^jlui |q| - w 2 
r(|q|- W ) 



(38) 
(39) 
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When q* < |q| < T, 



ReDiflql-w.lql.TJwdr + Ca 



Im£>i(|q|-a;,|q|,r)wC2 



where d = A^, C* 2 = ^ (i + 26- 1 ), C 3 = £L,C d = After direct calculations, we find 



64 i 



48tt ' 



J i K fa' 



7 2 oc 



^3 oc Vu/r. 

It is easy to check that 3> |/ 2 | ^ |^i|- Therefore, the temporal damping rate at T ^> lu has the form 

ImEL (w, T) oc vluT. 



(40) 
(41) 

(42) 

(43) 
(44) 

(45) 



This expression is surely not of the normal Fermi liquid type. It is, however, not the standard marginal Fermi liquid 
behavior. 

We next consider the transverse component of fermion damping rate. Similar to treatments presented above, the 
momentum |q| should also be divided into three sections, so that in the limit |k| = ,the damping rate IuiE,t(uj,T) 
has the form 



ImEx(w, T) 



1 

'2tt 



d|q||q|Im 



1 



D x (|q|-w.q,T) 



[nu(|q| -w) + ra F (|q| 



The approximate expression for £)j_(|q| — ui, q, T) were shown as follow. When < |q| < ^, 

ReZ)_i_ (|q| - |q| ,T) „ C.^^l! - C2 T(, Hq l)vf^M 



ImDx(|q|-w,|q|,T)wC3 



(w-|q|)Vw 2 -2w|q| 



When I < |q| < q*, 



ReDxdql-wJql.TjwCi 



T(\d\-u>Y , |q| y/2u |q| - ^ 



r 



ImD x (|q|-w,|ql,T) w-C a 



T(|q|- W )^2 W |q|- W 2 



(46) 

(47) 
(48) 

(49) 
(50) 



When g* < |q| < T, 



ReD ± (\q\-cj,\q\,T)^C 1 T 
ImDx(|q|-a;,|q|,T)«-C : 



TJ2uj 



(51) 
(52) 



whcrcC, 



4 247T 



After direct calculation, we can get 



Tl W 



In OC 



r 2 



^3 oc wloT 



(53) 

(54) 
(55) 
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It is clear that \I' 2 \ » » then we can conclude that, in the limit T » co 

T 2 

ImSx (w,T) oc — . (56) 

CO 

This expression is also divergent at zero energy u> — > 0. 

The above calculations were carried out in the zero-momentum limit |k| = 0. Using the on-shell approximation 
w = = |k| at finite temperature T 3> Uk, we find that ImSL (oj^,T) cx y/co-^T, but that ImEx is divergent 

in the infrared region. If we instead work in the zero-energy limit to = at T 3> |k|, then we find that Im^L (|k|, T) oc 
vlkjT and ImS T (|k|,T) oc ^/M?- 



B. Self-consistent computation of fermion damping rate 



The perturbative results for Dirac fermion damping rate always contain divergence, at both zero and finite tem- 
perature. It seems difficult to eliminate these divergences by including higher order corrections because they arise 
essentially from the singular gauge interaction. It is worth pointing out that such divergence exits in a wide range of 
physical problems. In the interacting electron gas, the electron self-energy diverges at zero energy if the bare, long- 
range Coulomb interaction is considered. In realistic metals, however, the bare Coulomb potential is always replaced 
by a short-ranged Yukawa potential after including the dynamical screening effect. The Debye screening ensures the 
infrared safety of the problem. In the present issue, however, the gauge interaction remains long-ranged even after 
taking the dynamical screening into account. Indeed, the gauge invariance ensures the masslessness of the U(l) gauge 
boson. Analogous divergence also appears in the perturbative computation for the damping rate of non-relativistic 
spinons due to scattering by gauge field [lj| . 

There is also an infrared divergence when computing pcrturbatively the self-energy of Dirac fermion in disordered 
potential (22j . Such divergence is essentially due to the linear spectrum of Dirac fermions. The most popular method 
to overcome this divergence is to invoke the so-called self-consistent Born approximation, which can give rise to a 
finite scattering rate [221 ] . It is natural to ask whether similar self-consistent approach can be used to eliminate the 
infrared divergence appearing in the present problem. To answer this question, we replace the internal free fermion 
propagator of Fig.l by the following retarded propagator 



G(w,k) 



— (j • k — ilm£(ui, T) ' 



(57) 



and then construct an integral equation for the damping rate ImS(w,T). If the perturbative expressions for the 
polarization functions are used, then only trivial result can be obtained after numerical computation. This is easy to 
understand by noting the fact that the infrared divergence originates from the singular gauge interaction. 

All the above computations are based on the polarization functions obtained from the free propagator of Dirac 
fermion. The feedback of fermion damping due to gauge field is completely neglected when calculating the dynamically 
screened gauge field and fermion self-energy function. However, if the fermion damping is really significant, as it 
should be in a non-Fermi liquid, its effect can not be simply neglected. It is conceivable to speculate that the 
divergences appearing in the Dirac fermion self-energy might be eliminated by incorporating the feedback effect of 
fermion damping. We now turn to such kind of self-consistent treatment of Dirac fermions and gauge bosons. Formally, 
the fermion damping rate satisfies the integral equation 



ImE(w,T) 



d 2 q 

xlm 



dcoi 



Im£(wi, T) 



7T ( Wl - |q|) 2 + (Im£( Wl ,T)) 2 



Di(wi-w,q,T) Dx(wi-w,q,T) 



[n B (u>i - lo) + n F (ui)] . 



(58) 



Using the fermion propagator (|57|) . the polarization functions ImDi and ImDi can be constructed as follows 



InxDr (e,|q|,T) 



d 2 k dcoi 



= 2N 



(2tt) 2 tt 
k • (k + q)Im 



Im 



lo\ — HwS(lo±, T) 
(wi -iIm£(a;i,T)) 2 - |k| 2 

Im 



Im 



{ui -ihnE(w 1 ,T)y - |k| 5 



Tm£(wi + e, T) 



(ui + e- ilmS(ui + e, T)) - |k + q| 2 
1 

(wi + e - !lmS(wi + e, T)) 2 - |k + q| : 
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x [np(uji) — np(u!i + i 



(59) 



ReD 1 (s,\q\,T) 

d 2 k duj 1 | 
2N / - — — - < Im 



(2tt) 2 tt 
-k • (k + q)Im 



UJ\ — iImS(o;i, T) 
(wi -iImS(o;i,r)) 2 - |k| 2 

1 



(wi - tImS ret (wi,T)) 2 - |k| 2 



x [n F (uJi) - n F (wi + e)] , 



Rc 



wi + e — ilmE (wi + e, T) 
(wi + e - ilmE (wi + e, T)) 2 - |k + q| 2 



Rc 



(wi + e - tImE rot (w! + e, T)) 2 - |k + q| 2 



(60) 



Im£ 2 (e,|q|,T) 

d 2 k duji 



= -AN 



(2tt)5 



Im 



wi — iImS(a;i, T) 
(wi -iImE(wi,T)) 2 - |k| 5 



x [n F (ui) - n_F(o;i + e)] 



Im 



(wi 



iImE(oji + e, T) 



iImE(wi +e,T))" 



(61) 



where 



Re£ 2 (e,|q|,T) 

. AT /" d 2 k du> lT 
-AN / - — ^ — -Im 

7 (2^) 2 TT 

x [njr(a;i) - n F {u-t + s)\ 



lji — iImE(wi, T) 
(wi -iImE(wi,T)) 2 - |k| 5 



Rc 



u!i + e — iImE(cji + e, T) 
e-ilm£(wi +e,T)) 2 - |k 



(62) 



Dx(e,|q| ) T) = Da(e > |q| J T) + T - i3 D 1 (e,|q|,r). 



(63) 



These equations appear to be very complicated and hard to solve. To get the fermion damping rate from these 
coupled inte gral equations, we find it is convenient to employ a simple dimensional analysis similar to that used by 
Vojta et al. [23j. After dividing all momenta, energy, and self-energy function by uj at T = and by T at finite T, 
it is found that ImE(a;) oc w and ImE(T) oc T respectively. Intuitively, this re-scaling procedure reflects the typical 
behavior of a marginal Fermi liquid. However, if we define 



ImE(w) 

ImE(T) 
f 



= A, 
= B, 



(64) 
(65) 



then numerical calculation find no convergent solutions for A and B. It turns out that, although formally the 
damping rate depends linearly on tu at T = and on T at uj = 0, a divergence appears in the regime where the 
energy of polarization functions vanishes. Specifically, this divergence emerges when lj\ — > lu in the fermion damping 
rate equation (|58|) . Since the functions D\ and D±_ appear in equation (|58[) as denominators, the fermion damping 
rate actually diverges as the energy of polarization functions vanishes. This qualitative analysis is confirmed by the 
numerical computations. 

In the above, we show that it is hard to get meaningful results of fermion damping rate within both perturbation 
theory and self-consistent treatment when the gauge field has no explicit Maxwell term. There is always some kind 
of divergence. All these divergences originate from the dynamically screened propagator of gauge field. The fermion 
damping rate seems not to be a well-defined quantity, at least under the approximations considered in the above. 



III. FERMION DAMPING RATE IN QED 3 WITH MAXWELL TERM 



In this section, we study the QED3 theory with explicit Maxwell term. According to the standard framework of 
rclativistic quantum field theory, it is natural to keep such kinetic term in the Lagrangian. The QED3 with a Maxwell 
term itself is very interesting and has been studied extensively in the past twenty years (for a review, see 24j ) . It also 
has direct applications in condensed matter physics. In the context of underdoped high temperature superconductors, 
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an effective QED3 theory was derived to model the unusual physics after carefully considering the phase fluctuations 
[uj[ll|]. There is a Maxwell term in this kind of QED 3 theory. In the following, we will include the Maxwell term for 
gauge field in the Lagrangian of QED3 and re-calculate the fermion damping rate. 
In the presence of Maxwell term, the propagators for gauge field are 



D 00 (n m ,q) 



Di(fi m ,q) 



QiQj 



1 



nl + D ± (fi m ,q) 



(66) 
(67) 



After pcrturbativc computations, we find the following transverse damping rate 
ImSx (w, k) : 



AN 



d z q- 



(w-|k + q[ 



iV 2 (( W -|k + q|) 2 -|q| 2 



256 ^-|k + q|) z -|q 
xe(w-|k + q|-|q|). 

Since 9 (w — |k + q| — |q|) = when \ui\ < |k|, we still have ImEx (uj < |k| , k) = 0. In the limit |k| = 0, we get 



In the low-energy regime, u) — > 0, it reduces to 



1 ,16w, N 

— arctanl ) H 

4tt N 64ttu; 64nuj 16w 



N N 16w 
arctan( — ). 



ImEx(w) 



8H 

3A/V 



(68) 



(69) 



(70) 



which coincides with the standard behavior of marginal Fermi liquid. The general ImSx (w, k) can only be calculated 
by numerical methods, with results shown in Fig. 01 





(a) (b) 
FIG. 4: Transverse contribution to damping rate for Dirac fermion in the presence of Maxwell term. 
On the other hand, the longitudinal damping rate is found to be 

AN 



ImS L (w, k) = =- / d'q 



( W -|k + q |) 2 -|q| 2 



256 (( W -|k + q|) 2 -|q| 2 ) |q| 2 + iV 2 |q| 
It is easy to prove that ImEL (w < |k| , k) = 0. In the limit |k| = 0, 

87V n 



-0( w -|k + q|-|q[ 



Im£ L (u) 



t Jo Iql 



- ^ |q| 



256 (w 2 -2u;|q|)+jV 2 



(71) 



(72) 
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This is still divergent in the infrared region. The general IhiSl (w,k) ban be written as 

4N f +oc 

Im£i>,k) = T d|q|F 2 (|q|) 

n Jo 



where 



F a (|q|)= / d v 
Jo 

When oj > Iql, we have 



( W -|k + q|) 2 -|q| 2 



256 (( W -|k + q|) 2 -|q| 2 ) |q|+iV 2 |q| 



i( w -|k + q|-|q[ 



(73) 



(74) 



lim F 2 (|q|) 

|q|->0 



1 



2tt (w- |k|) 
256 (u - |k|) 2 + N 2 



(75) 



Thus we can easily see that Im£L(w,k) is infrared divergent when |w| > |k|. 

Due to this divergence, the longitudinal fermion damping rate is still ill-defined at zero temperature T = 0, even in 
the presence of the Maxwell term. 

At finite temperature, after tedious calculations we find that the Dirac fermion damping rate has the following 
features: If we calculate in the limit |k| = at T w, then Im£L (w, T) oc \fuff and ImEx (w, T) oc T/uj. If we work 
in the on-shell approximation u = = |k| at T lo^, then we find that Im^L (wk, T) oc y/uj^T and that Im£x (^k, T) 
is divergent in the infrared region. In the zero-energy limit oj = with T 3> |k|, we have Im^L (|k|,T) oc y/\k\ T and 
Im£ T (|k|,T)oc y/\k\T. 

All the above calculations are obtained in the Coulomb gauge, which separates the longitudinal and transverse 
components of the gauge field completely. The same calculations can be done similarly by choosing another gauge. 
After straightforward computation (25| . we find that the fermion damping rate is still divergent in a general gauge 
when obtained at the pcrturbative level, no matter the Maxwell term of gauge field is present or not. This implies 
that the existence of divergence in perturbative expansion is a universal feature of QED3, rather than just a gauge 
artifact. 

The marginal Fermi liquid behavior of fermion damping rate was claimed previously by Franz and Tesanovic without 
providing computational details [loj . However, the detailed calculations show that the fermion damping rate is an 
ill-defined quantity at the perturbative level because divergence appears at both the T = and T>w limits. 

We next turn to the self-consistent treatment of fermion damping rate in the presence of Maxwell term. To compare 
with the results presented above, we also choose to work in the Coulomb gauge. As in the last section, we include 
both the real and imaginary parts of the vacuum polarizations when writing the integral equation for the fermion 
damping rate 



Im£(w,T) 



1 


p + OC 




Jo 


•{ 


-Im 



duii 



1 



Im£(oJi, T) 



(wi- |q|) 2 + (Im£( Wl ,T)) 2 



Dx (wi-wJql.T) 



Im 



(wi-wr + D^ (wi-wJql.T) 



x [ns(wi -uS) + n F (uji)] 



(76) 



where the functions Di and D± are given by equations (|59H63[) in the last section. The kinetic term of gauge field 
eliminated the divergence brought by the polarization functions. At the same time, the energy uj can no longer be 
scaled out due to the kinetic gauge term. Thus the fermion damping rate is not expected to display the marginal 
Fermi liquid behavior. After numerically solving these coupled integral equations, we find that the total fermion 
damping rate is 



ImE(w) 



5( a ; 



(77) 

The computation of fermion 



at zero temperature T = 0. The numerical results at T — arc presented in Fig. 
damping rate at finite temperature follows the same procedure as presented above. The numerical computations find 
the following fermion damping rate 



Im£(T) oc VT 



(78) 
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in the limit T ui. The numerical results for this limit are presented in Fig. 5(b) As a summary, the damping rate 
of masslcss Dirac fermions is 



ImS(cj, T) oc max 



T) 



(79) 



due to scattering by the U(l) gauge field in two spatial dimensions. This is certainly a non-Fermi liquid like behavior 
since the fermion damping rate would be ImS(w, T) cx max(w 2 , T 2 ) (or with higher powers) in a normal Fermi liquid. 






0.004 0.006 
Tla 



(a) 



(b) 



FIG. 5: (a) Damping rate of massless Dirac fermions at zero temperature T — 0. The fitting curve is ImE (u) /a = 
0.2170(a;/a) ' 5051 . Here the energy scale a is defined as a = Ne 2 /8 (we have let e = 1 in this paper), which was first in- 
troduced in Ref. [l7j |. (b) Damping rate of massless Dirac fermions on Fermi level at finite temperature. The fitting curve is 
ImE (T) /a = 1.4220(r/a) ' 5014 . 



This treatment is essentially the analogy of Eliashberg theory of superconductors with strong electron-phonon 
interaction. The propagators of fermions and intermediate bosons are self-consistcntly coupled while the vertex 
corrections are simply ignored. In the standard Eliashberg theory of electron-phonon system, the vertex corrections 
can be safely neglected since they are suppressed by the small parameter m/M, where m and M are the electron 
mass and nuclei mass respectively. This is nothing but the Migdal theorem [261 ]. In the present problem, there is no 
similar mass scales since both Dirac fermions and gauge bosons are massless. However, we do have a small expansion 
parameter 1/N. The vertex corrections are suppressed by the factor 1/N. In the large N limit, our ignorance of the 
vertex corrections is valid (Polchinksi used the same argument when studying the fermion self-energy due to gauge 
interaction, see his paper in Ref. [HI). The validity of our results for small N, such as N = 2, is not clear currently. 

Unfortunately, at present we are unable to get a well-defined result for the momentum dependence of fermion damp- 
ing rate from the same self-consistent treatment. The reason is that the coupled equations become very complicated 
after including the momentum dependence. The numerical results are much less reliable as those at zero momentum. 



IV. MARGINAL FERMI LIQUID BEHAVIOR AND COULOMB INTERACTION 



The marginal Fermi liquid was proposed by Varma et al. at 1989 to understand some of the highly unusual 
experimental facts in the normal state of high temperature superconductors on phenomenological grounds [21| . Since 
then, a lot of efforts have been devoted to deriving this phenomenological theory from certain microscopic models 
0, E3, EE EH • Gauge theory has long been considered as one possible candidate [l(| El, Ell ■ Within the ordinary 
perturbative theory, the transverse damping rate of Dirac fermions has a linear dependence on energy at T = 0, which 
is the standard behavior of marginal Fermi liquid. However, this result can not be trusted because of the appearance 
of divergence in the longitudinal component at T = and in the transverse component at finite temperature. More 
careful theoretical and numerical computations show that the Dirac fermions exhibit non-Fermi, rather than marginal 
Fermi, liquid behavior. 

However, it is possible to find signature of marginal Fermi liquid within some models those are in form analogous 
to the U(l) gauge field theory. For instance, we now consider the long-range Coulomb interaction between Dirac 
fermions, which remains unscreened because of the vanishing density of states at the Fermi level. In some sense, the 
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Coulomb potential may be regarded as the temporal component of a U(l) gauge field. There is, however, a subtle 
difference, which leads to important consequence. The dynamically screened Coulomb interaction can be formally 
written as 

D c (in m ,q) = -, (80) 

if +Di(ifi mj q) 

where the |q| term comes from the bare, instantaneous Coulomb potential and the dimensionless parameter is defined 
as A — — with eo being the dielectric constant. At zero temperature, the damping rate of Dirac fermions has the 
form 



4N r , v ( w - l k + q|) 2 - |q| 2 

ImS c ( W ,k) = -— d\ V ' q w _ k + q|-|q|) (81) 

" J 256 (( W -|k + q|) 2 -|q| 2 )^ + 7V2|q| 2 

Here, we also have ImEc (w < |k| , k) = 0. In the limit |k| = 0, we get 

Im£ c (u) = -^C{\)u} (82) 

with 



C{\)= [ dx |VL^ . ( 83 ) 

This damping rate is free of divergence and can be identified as the behavior of a marginal Fermi liquid. It agrees 
with the results obtained in the context of graphene [27M29j . 

The damping rate ImSc (w, k) at general energy-momentum is calculated by numerical methods, with results shown 
in Fig. [6] 

We can also study the self-consistent equations of damping rate and polarization function. As in the case of gauge 
field, we consider only the zero-momentum limit. Now the equation for fcrmion damping rate is 



1 f + °° , , f + °° haS(ui,T) 

ImS (w,T) = -— =■ / d q / dui )-J±-i Im 

2tt 2 7 IMI i-oo {uJi - |q|) 2 + (ImS(wi,T)) 



1 



q| 



, Di(wi-w 1 q J T)_ 

x[n B {wi-u) + n F {ui)], (84) 

where D\{e, q,T) was given by (|59|) and (|60|) . Due to the special form of the bare interaction function |q| _1 , the 
energy w can be completely scaled out in the whole set of integral equations. This is not possible for |q| -2 , as discussed 
in Sec. 3. Moreover, such bare term ensures the absence of divergence. In the absence of any bare term, the re-scaling 
can still be performed, but the equation is divergent, as discussed in Sec. 2. Therefore, the Coulomb interaction is 
very special and turns out to be a good candidate for producing marginal Fermi liquid behavior. 

With the help of scaling analysis, it is easy to show that ImS(w) oc 10 in the limit u ^> T and ImS(T) oc T in the 
limit uj <C T, respectively. As before, we define Im ^(") = A and Im ^( T ) = ^ then their dependence on A can be 
obtained by numerical computation. The results are presented in Fig. 7(a) and Fig. |7(b)[ 



The linear dependence of fcrmion damping rate on energy /temperature may appear in other counterpart of the 
U(l) gauge interaction. It was discovered in the model of nodal quasiparticles coupled to critical fluctuation of some 
superconducting order parameter [23l. [301] . This model, albeit having different physical contents, share one common 
feature with Coulomb interaction: the coupling of massless Dirac fermions to some singular boson mode. This is 
in form analogous to U(l) gauge interaction. However, the calculation of fermion damping rate caused by gauge 
interaction meets with divergences. According to our theoretical and numerical computations, it seems that the 
kinetic term for gauge field has to be explicitly included in order to get meaningful results. 

Aji and Varma constructed an interesting dissipativc quantum 2D XY model and showed that it produces marginal 
Fermi liquid behaviors [3l| . Recently, behavior of marginal Fermi liquid type was also found in a d-dimensional field 
theory with the help of AdS/CFT correspondence (32j . 

We also note that a linear-in-T quasiparticle damping rate was already pointed out in two early papers [3, • I n 
both of these works, the linear-in-T behavior is attributed to the scattering of quasiparticiles by an emergent U(l) 
gauge field. However, such damping rate is that of spinless holons, rather than fermions. Therefore, the interacting 
system considered in these papers, though very interesting, can not be identified as a marginal "Fermi" liquid. 
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(c) (d) 

FIG. 6: Damping rate for Dirac fermion due to the Coulomb interaction. A" 1 are 0.001, 0.01, 0.1, and 1 in (a), (b), (c), (d) 
respectively. 

V. SUMMARY AND DISCUSSION 

In summary, we present a detailed calculation of the damping rate of massless Dirac fermions due to U(l) gauge 
field in QED3. When the theory contains no Maxwell term for gauge field, the fermion damping rate is found to 
diverge at both zero and finite temperatures within perturbation theory. There is still divergence in the self-consistent 
equations for fermion damping rate and gauge boson propagator. Once the Maxwell term for gauge field is included 
into the self-consistent equations, the fermion damping rate is well-defined and display non-Fermi liquid behaviors at 
both zero and finite temperatures. 

From the first sight, the existence of divergence in fermion damping rate in QED3 theory without Maxwell term 
might restrict its validity in understanding high temperature superconductors. However, in reality this is not as severe 
as it looks. Such effective theory applies to both the half-filling state and the underdoped region of high temperature 
superconductors. In the half-filling state, the massless Dirac fermions undergo a pairing instability towards the chiral 
symmetry breaking phase. Now the fermionic excitations are suppressed by the dynamically generated mass gap and 
thus it is usually not necessary to study the damping rate. In the underdoped region, the gauge field couples not only 
to massless Dirac fermions, but also to an additional boson field which describes the motion of charged holons. The 
interaction between gauge field and holons contributes a vacuum polarization function to the gauge boson propagator, 
which might be able to eliminate the divergence appearing in the fermion damping rate. To address this issue, it is 
essential to carefully study the whole interacting system, especially the holon-gauge coupling. 

At present, we are unable to eliminate the divergence appearing in the momentum dependence of Dirac fermion 
damping rate. Such divergence might also be cured by the self-consistent (Eliashberg) treatment, but the coupled 
equations become much more complicated than the zero-momentum limit and are hence hard to be solved numeri- 
cally. We will study this problem further, either by improving numerical methods or by developing novel theoretical 
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approaches. 

When the U(l) gauge theory is used to describe the anomalous properties of high temperature superconductors, the 
Dirac fermions are usually not the physical electrons. They might be fermionic spinons 0, [B43, fermionic holons 
or topological fermions pH EH , depending on the physical motivations and the procedures of deriving the effective 
field theory. Therefore, the damping rate studied in this work could be directly compared with experiments only after 
including the additional degrees of freedom. However, since the U(l) gauge interaction of massless Dirac fermions 
appears naturally in a number of correlated electron systems, we believe it is interesting to carefully study the damping 
rate and other physical quantities of the "unphysical" Dirac fermions. 
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